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Abstract: Relevant deformations of gravity present an exciting window of opportunity 
to probe the rigidity of gravity on cosmological scales. For a single-graviton theory, the 
leading relevant deformation constitutes a graviton mass term. In this paper, we investigate 
the classical and quantum stability of massive cosmological gravitons on generic Friedman 
backgrounds. For a Universe expanding towards a de Sitter epoch, we find that massive 
cosmological gravitons are self-protected against unitarity violations by a strong coupling 
phenomenon. 
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1. Introduction 

Technical naturalness is arguably one of the most promising pathfinders to physics beyond 
the standard model of particle interactions and gravity, as well. It offers many exciting 
windows of opportunity related to renormalisable standard model operators that share an 
enhanced sensitivity to the scale of new physics. 

Among these, the vacuum energy density is standing out in various ways. Being the 
unique operator with quartic sensitivity to the ultraviolet scale, it also represents the most 
relevant term in the Einstein-Hilbert action. The basic observation is that the vacuum 
energy density is technically unnatural within the standard model already at energy scales 
set by the lightest measured particle masses within its spectrum. In other words, the 
technical naturalness facet of this challenge is not solely tied to the quantum gravity scale, 
unless there is an ultraviolet-infrared conspiracy operative in the vacuum sector that also 
respects the many high-precision successes of the standard model at lower energies. 

Although collider experiments cannot measure the vacuum energy density, the chal- 
lenge it poses becomes serious once the standard model of particle physics is coupled to 
gravity. Since gravity, being the most democratic field theory, couples to energy-momentum 
in a universal manner irrespective of its sources' nature, the vacuum curves spacetime and 
affects the Universe's expansion history. As a consequence, and opening up yet another 
window of opportunity, it seems attractive to reconsider gravity in the deep infrared and 
its consistent deformations. 

The hunt for a fundamental completion of gravity is characterised by incorporating new 
degrees of freedom in the ultraviolet with the infrared region kept untouched (in a relevant 
sense), its role reduced to providing the classical benchmark tests. However, this precludes 
the opening up of additional gravitational degrees of freedom that could be relevant in the 
deep infrared on scales that are only poorly constrained by state-of-the-art cosmological 
observations. 
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Historically, this question was discouraged by a mighty no-go theorem stating the im- 
possibility to embed gravity in a QCD-like theory with a self-interacting graviton multiplet 
under the spell of Yang-Mills. More precisely, assuming locality, Poincare invariance, and a 
free-field limit consisting of massless gravitons, the only consistent deformations involving 
a multiplet of gravitons are such that the deformed gauge algebra is just a direct sum of 
independent diffeomorphism algebras Q. 

Interestingly, by relaxing one of its conditions — allowing for new relevant degrees of 
freedom — this no-go theorem gave rise to a potent tool for studying consistent defor- 
mations of Einstein's gravity within the effective field theory framework. For a nonzero 
deformation parameter, the candidate deformations are characterised by a tremendous 
reduction of symmetries: 

diff(Mi) ® diff(M 2 ) ® > diag(diff(Mi) ® diff(M 2 ) »•••) , 

with Mj denoting (not necessarily) different spacetime (sub-)manifolds. The resulting 
gauge symmetry restricts the deformation, i.e. the coupling of different geometries, to be 
solely constructed from invariants of the reduced symmetry ||] . 

Under the umbrella of this framework, many proposals for relevant deformations of 
gravity that have been suggested in the last decade become cousins. Even more promising, 
the unique ghost-free theory for a massive spin-two field propagating all degrees of freedom 
H fits under this umbrella. In the latter case, two copies of Minkowski spacetime are 
considered, one perfect the other perturbed, and the most relevant deformation becomes a 
spin-two mass term 

5 deform = ~ J d 4 x V • (1-1) 

The mass matrix M depends only on the background geometry and is constant in this 
case. As mentioned earlier, this matrix is uniquely determined by unitarity arguments and 
requires tuning. 

There is a rather straightforward nonlinear completion of the leading infrared defor- 
mation Q, 

S de{om = ~ ! d*xV=FiM a ^( gi ) H a pH„ u ( gi ,g 2 ) , 

1 J Mi 

Ha/3 = 9la/3 — E^Ep g 2 ftu , (1-2) 

where E denotes the pullback from M 2 to M\. Since the spacetimes (Mi, g\) and (M 2 ,g 2 ) 
need not be diffeomorphic to each other (even not in the perturbative sense), H{x) is in 



general not a fluctuation on M\. In fact, the deformation (1.2) represents a mass term for 
a graviton 1 only if (M 2 ,g 2 ) is a copy of (Mi,gi) at the background level. 

Consider the case when the most relevant deformation is a mass term. Let us assume 
for simplicity, that some sort of Higgs mechanism is responsible for the graviton mass gener- 
ation. For instance, in the Fierz-Pauli setup, a massless graviton, which has two transverse 



1 In abuse of notation and logic, we refer to metric perturbations as gravitons irrespective of the back- 
ground geometry. 
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polarisation states, combines with a Goldstone vector to a massive spin-2 field, which has 
six polarisation states in general (but only five on Minkowski spacetime). The Goldstone 
vector carries three transverse polarisation states and one longitudinal polarisation state. 
When the massive graviton is at rest, its six polarisation states are completely equivalent. 
However, if it is moving, the longitudinal polarisation becomes increasingly parallel to the 
graviton's momentum. As a consequence, at high energies, a massive graviton might look 
like the longitudinal polarisation state carried by the Goldstone vector 0]. This state- 
ment is known as the Goldstone equivalence theorem and rests on the underlying gauge 
invar iance. 

The longitudinal polarisation state does not receive kinetic support from the Einstein- 
Hilbert term, which is precisely why it is in the focus of unitarity requirements. As a 
matter of fact, many distinguished features of massive gravity, like, for example the van 
Dam-Veltman-Zakharov discontinuity Q j|] , the Vainshtein radius 0], or the structure 
of the Fierz-Pauli mass term on Minkowski spacetime ||, are captured by the longitu- 
dinal polarisation state of the Goldstone vector. Consistency invesitgations of massive 
gravity mostly focused on ghostly excitations at some finite perturbation level based on a 
Minkowski ground state, with the earliest exception being Higuchi's unitarity analysis on 
de Sitter spacetime ||. 

Higuchi found a consistency relation between the deformation parameter, the graviton 
mass, and the curvature scale of de Sitter, set by the cosmological constant: In order to 
avoid negative norm states m 2 > H 2 , where H denotes the Hubble constant. This bound 
is of great interest, since de Sitter geometry is unique in the sense that is does not require 
any source specification. However, from a field theoretical point of view this makes the 
setup special, because the background reference scale is constant here. 

In this paper, we generalise Higuchi's bound from de Sitter to general Friedman cos- 
mologies by employing the Goldstone equivalence theorem outlined above. We find a 
competition between classical stability, the requirement that perturbations respect the 
background, and quantum stability, the requirement that the spectrum does not contain 
negative norm states. For the special case of a de Sitter background, both criteria coin- 
cide and give rise to the unitarity bound quoted above. The situation is richer for generic 
Friedman cosmologies. There, the very nature of either bound is more intriguing, since it 
involves a time dependent curvature scale that is monotonously increasing or decreasing in 
the past, depending on the specific sources that drive the background expansion. Generi- 
cally, none of the bounds can be satisfied on the entire spacetime manifold. This, however, 
does not imply that the theory is invalidated. Indeed, it turns out that massive gravitons 
in generic Friedman universes are protected against unitarity violations. More precisely, 
it is not clear at all whether a unitarity bound really exists, because before entering the 
would-be unitarity violating spacetime region, the theory becomes strongly coupled 2 . 

We reach the following verdict: Phenomenological constraints require to choose the 
initial hypersurface close to the present hypersurface. In the case of a radiation or matter 
dominated Friedman universe, the evolution towards future hypersurfaces is guaranteed to 

2 The connection between stability and strong coupling is clarified at the end of section H. 
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be healthy (for consistent initial conditions) by the strictly monotonic background expan- 
sion. In the most interesting case of a matter-cosmological constant mixture, the future 
evolution will be sound provided the mass is large enough and consistent initial condi- 
tion have been imposed. In all of the above cases, evolving backwards in time, massive 
cosmological gravitons will soon enter a strong coupling regime that demands a nonlinear 
completion of the theory. In other words, the fact that we have a sound theory on all past 
hypersurfaces is nontrivial and ensured by a strong coupling phenomenon that is confined 
to these spacetime regions. This is the advertised self-protection mechanism. 

Last but not least, we show that all conclusions hold for a generic Friedman source. 

2. Goldstone— Stiickelberg analysis 

In this section we derive the classical and quantum stability requirements for massive 
cosmological gravitons. 

Consider two copies of a generic (background) spacetime, (Mb, 7) and (Mb, g) with 
7 = g + h, where h denotes the metric perturbation obeying \g(t)\ 3> |/i(i,x)|. In this case, 
H = h is a perturbation (where, for simplicity, we have chosen the same coordinate system 
on both manifolds) under the spell of diff(MB) for vanishing deformation parameter. 

Turning on the most relevant deformation (|1.2| ) , the gauge symmetry is deformed to the 
diagonal subgroup of diff(MB) <8> diff(MB). In other words, the deformation removes the 
freedom to gauge h relative to the background geometry. The massive graviton, however, 
still carries six degrees of freedom, due to the second Bianchi identity 

= V u h . (2.1) 

In general, this is not a gauge choice. For instance, in the undeformed theory for massless 
gravitons on Minkowski spacetime the constraint fl2.1|) is not a legitimate gauge, because 
the corresponding gauge shifts become singular for this choice (as a testimony of the van 
Dam-Veltman-Zakharov discontinuity on this background). 

As is well known, in view of this explicit symmetry deformation, there are two equiv- 
alent state descriptions. In the first case, the metric perturbation is split according to 

V = h U + V (mK) , (2-2) 

where h is covariantly conserved and carries two transverse degrees of freedom, while V is 
unconstrained and carries four degrees of freedom. The latter can be decomposed further, 
V = V 1 - + d^> . Here V • V 1 - = and Vl/ carries one degree of freedom. Using this state 
description, the theory can be considered as a gauge fixed theory. 

The second and equivalent state description, called the Goldstone-Stiickelberg com- 
pletion, is based on adding four degrees of freedom, carried by a vector field ir in order to 
restore the original gauge symmetry, diff(MB) <8> diff(MB). In this case, the completion is 
given by 

Hfj, u = + V( M 7Tj,), (2.3) 
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where H has ten degrees of freedom, of which six are carried by h and four by the Goldstone- 
Stiickelberg vector it, which can be further decomposed as tt = ir 1 - + dcf). Here V • tt 1 - = 
and (j> carries one scalar degree of freedom. 

As mentioned above, the crucial point of this construction is the restored gauge sym- 
metry that allows to shift h and tt relative to the background such that the Goldstone- 
Stiickelberg completed H itself is rendered gauge invariant. It is clear that four degrees of 
freedom represent gauge redundancies, leaving us with six physical degrees of freedom. 

We choose to work with the Goldstone-Stiickelberg completed state, for which the 
leading relevant deformation is exactly the celebrated Fierz-Pauli mass term, 

,2 



Mr 



Mafinvig) = ga^gpu - 9ap9^ ■ (2.4) 

The mode of the metric fluctuation corresponding to the Goldstone-Stiickelberg scalar 
</> dominates scattering processes at high momenta. This is tantamount to the Goldstone 
boson equivalence theorem and most easily understood from the observation that this mode 
enters processes with at least two derivatives, W^d u (p, which, therefore, grows fastest in the 
high momentum limit. This is precisely the regime for which we are interested in studying 
the stability of the deformed theory. 



The field 4> mixes with the metric perturbation h through the mass term (2.4), 



2m 2 (Dh - VVV) = 2m 2 (r<$ (g)h^ - i?« (g, hj) <t> , (2.5) 



where we have integrated by parts and introduced the Ricci tensor R^ = R evaluated on 
the background configuration g, as well as the Ricci scalar R^ expanded to first order in 
h. 

In order to eliminate the kinetic mixing term R^(g,h)4>, we carry out a conformal 
transformation 

9fip = ^9iiv = (1 + w) 2 ^, (2.6) 

which, at the linear level, is equivalent to = hn V + 2wq^ u . The Einstein-Hilbert term 
transforms as 

^gtt 2 R = y/^§ (r - 6JT 2 g ab d a ttd b n) . (2.7) 



In order to eliminate the mixing between h^ u and <f> we must choose Q 2 = 1 — 2m?(f) or, 
equivalently (since is a first order quantity in the expansion of H^ u ), lo = —m 2 (j). 

The conformal transformation ( p. 7] ) contributes a standard kinetic term for the Goldstone- 
Stiickelberg scalar, while the massive deformation fl2.4p gives rise to a non-standard kinetic 
contribution with the metric field replaced by the background Ricci tensor, 

2m 2 ((D^) 2 - V^V^VVV) = 2m 2 R llv d' i (t)d u (t) . (2.8) 

The action for <p is given by 



S = J d 4 x^g (A<p 2 + B^(di(l))(d j( f>) + j>D%<fy , (2. 
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with A = 2m 2 (-3m 2 g 00 + R 00 ), = 2m 2 (-3m 2 g i i + R ij ), and D l = Am 2 (-3m 2 g i0 + R i0 ) . 
In general, these coefficients are spacetime dependent functions. Note that in ( |2.9| ) we 
have not displayed any potential terms (self-couplings) such as R(g)hcf), since quantum 
stability refers to the free evolution, and classical stability relies on the kinetic terms at 
high momenta. The corresponding Hamilton density reads 

n = -iU - x/^WXc^) , (2.10) 



A^gA 

where tt = 5£/<f). The Hamiltonian is unbounded from below for A < or B lJ positive 
definite. 

In the case of a generic Friedman background geometry the action for the Goldstone- 
Stiickelberg scalar reduces to 

A(t)<j> 2 + B(t) (v<f>/aY) , (2.11) 



S 



d 4 ^ 



-g 



where A(t) = 6m 2 (m 2 — H — H 2 ), B(t) = -2m 2 (3m 2 - H - 3H 2 ), H = H(t) denotes 
the Hubble parameter, and o = a(t) the scale factor. These coefficients depend on the 
energy-momentum source curving the Friedman background and, in particular, can change 
signs during the background evolution. Classical stability requires A > and B < 0. 

In order to study quantum stability of the Goldstone-Stiickelberg scalar on a generic 
Friedman background, we need to transform ( 2.1l| ) into normal form. This requires the 
following field redefinition: </> -»• fcj) with / = -Cf/(2A), where C = A- 3HB. The 
Lagrangian transforms as 

f~ 2 C = A ft - C # + B (V<p/a) 2 + C 2 /(4A) (f? . (2.12) 



The coefficient A controls the sign between (j) and / 2 n = 2A(p — C(j) and has therefore 



an important impact on the quantum stability of (2.11). This can be worked out along the 
canonical quantisation prescription. As usual, we postulate the equal time commutation 
relations 

[0(i,x),7r(t',x')] t=t/ =i5 {3) (x-x') , (2.13) 
and decompose the field (j)(x) into modes U(t, k) = u(t, k) exp (ik • x), where u(t, k) satisfies 

2 



Au 



a 



d C\ C 
+ 1 di + 2A 2A 







(2.14) 



The fact that the collection of U(t, k) represents a complete orthonormal set of solutions 
(with respect to a spatial hypersurface-independent scalar product) results in a simple 
condition on the Wronskian of the solutions, 



(u*ii - uii*) (t,k) = 1 . 

As a consequence, the Goldstone-Stiickelberg scalar may be expanded as 
d 3 k 1 



(2.15) 



(2^)3/2 ^2\A\f 



u(t, k) exp (ik • x)a(k) + u*(t, k) exp (— ik • x)a^(k) ) . 



(2.16) 
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Inserting this expansion and the corresponding one for tt into the canonical commutation 
relations ( |2.13| ) yields 

a{k),a\k')\ =sign(A) 5 (3) (k-k'). (2.17) 

The construction of a vacuum state and Fock space can now proceed as usual. However, 
whenever A < 0, the construction results in negative norm states, which violate unitarity. 

For an arbitrary spacetime, this quantisation procedure bears conceptual challenges, 
because there may be no Killing vectors at all to define positive frequency modes. The 
situation is simpler for a Friedman background since it accommodates a restricted set of 
isometries, i.e. invariance under spatial rotations. Then, together with the corresponding 
Killing vectors there exist associated (natural) coordinates. Of course, coordinate systems 
are physically irrelevant — a fact that renders the particle concept somewhat arbitrary on 
curved spacetimes. However, this concerns the interpretation of the theory. The unitarity 
requirement A < is a coordinate independent statement, since Bogolubov transformations 
are unitary transformations. 

The most important result of this section is that stability considerations led us to 
require A > and B < in order to have a bounded Hamiltonian and, in addition, A > 
to have a sound probabilistic interpretation. Hence, provided the theory respects unitarity, 
B < represents the classical stability bound. 



This is also clear from the equation of motion (2.14). Indeed, for A > 0, to have a 
stable solution, requires B < 0. Further, in order to have a damped solution at late (early) 
times demands C > 0(C < 0). Whenever this condition is violated, the background will 
be destabilised. 

In order to investigate the signs of the coefficients A,B,C in a spatially flat Friedman 
Universe, consider the Friedman equations 

3H 2 = 8irGp + A, (2.18) 
3 (H + H 2 ) = -A-kG {p + 3p) + A, (2.19) 



where p is the density, p the pressure, and A the cosmological constant. An expanding 
universe is characterised by H > and p < 0, so H < 0. Furthermore, H = —3HH{l+c s 2 ), 
where c s denotes the isentropic sound speed. As a consequence, H > 0, for an equation of 
state governing an arbitrary mixture of matter and radiation. 

In the absence of a cosmological constant, H and all its time derivatives vanish at 
late times. The late time asymptotics of the coefficients are therefore given solely by the 



gravitons mass, A = Qm 4 = —B > and C = 0. Consequently, the modes ( 2.14 ) are 
well-behaved. In a de Sitter universe, A = 6m 2 (m 2 — H 2 ) = —B and C cc A. The absence 
of negative norm states requires m 2 > H 2 , which is nothing else but the famous Higuchi 
bound. We see that in this case, classical and quantum stability collapses to a single 
criterion between the graviton mass and the constant de Sitter curvature scale. 

For a generic Friedman universe, at any moment in time, A > —B and C = 6m 2 H(3m 2 — 
H/H — 3H — 3H 2 ) < 3HA. As a consequence, evolving the modes ( |2.14| ) backward in time, 
B and C change signs before A does (i.e. at later cosmological times). In the case of C this 
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sign change stabilises the modes since it leads to mode damping. This stabilisation, how- 
ever, is only marginal and nullified by the change of sign of B which triggers an exponential 
instability that dominates the large proper momentum regime. 

As a consequence, evolving the modes backwards in cosmological times, the system 
enters first a strong coupling regime (at later cosmological time) before it would violate 
unitarity (at earlier cosmological times). In this sense, the strong coupling regime protects 
massive cosmological gravitons from unitarity violations. In other words, whenever massive 
gravitons in an expanding Friedman universe experience unitarity violations, then, for sure, 
they are already in a strong coupling regime that demands a nonlinear completion, but not 
vice versa. In this respect, de Sitter is a borderline geometry since both inconsistencies 
coincide. 

It is important to appreciate that once a generic mode enters the classical instability 
region B > 0, it does not destabilise the background instantaneously. Instead, the charac- 
teristic time scale for this to happen is T(k) oc a/k. As a consequence, modes with arbitrary 
large proper momentum become strongly coupled without further delay once B > 0. More 
precisely, for arbitrary initial conditions, there will always be a critical proper momentum 
k* such that all modes with momenta k > k* enter the nonlinear regime before they would 
violate the unitarity bound. 

Strictly speaking, the self-protection mechanism is confined to and efficient for B < 0, 
that is, before the fluctuations enter the unstable regime, which is characterised by the 
background coefficient B becoming positive, they hit a strong coupling regime. Indeed, for 
B < and B — > 0, the proper kinetic energy density oc B(S7<j)/a) 2 becomes subdominant on 
any scale as compared to the potential energy. Hence, the exponential instability discussed 
above, characterized by B > 0, is a testimony of this strong coupling regime and can 
legitimately be used to identify it. 

2.1 Self-protection under matter impact 

In the following we argue that both the stability and unitarity behaviour of the theory 
are generically not altered by the specific choice of a matter action Setter r&% > g^v] , where 
denotes the collection of matter degrees of freedom. This is remarkable because the 
matter action explicitly depends on the metric field and, in particular, on the Goldstone- 
Stiickelberg scalar <fi. This coupling could, in principle, change the dynamics of <j) in a 
relevant way such that the self-protection mechanism will be overridden. As we will show 
below, this is not the case. 

Since the matter action is invariant under general coordinate transformations, the 
Goldstone-Stiickelberg scalar enters the matter sector only via the conformal transforma- 
tion fl2.6|) . Let us first expand the matter action to second order in the fluctuations, 

Matter [*< + 6*U9 + h] D j d 4 X d 4 y h^(x) U a ^y)^J— + T ^ (x) . 

(2.20) 

After the conformal transformation fl2.6|) , h = h + 2m 2 (f)g, the first term on the right-hand 
side of (2.20) contributes only self-interaction terms of the form </> 2 and couplings h(j>. The 



-8- 



coefficient of this potential term can be estimated to be of order H 2 . Hence, the potential 
term is subdominant as compared to the kinetic term at high proper momenta. 



The conformal transformation of the second term on the right-hand side of (|2.2C| ) 
gives rise to a coupling between the Goldstone-Stiickelberg scalar and the matter degrees 
of freedom, 

/BT 
d ' x ^WM 3 ^ 1 ' (2,21) 



where we have again neglected potential terms of the form 4>5^i. A coupling like ( 2.21 ) 



will not modify the momentum field conjugated to </>, and, hence, the unitarity bound is 



robust against its inclusion. Suppose now we integrate ( |2.2l|) by parts, thereby producing 
a time derivative acting on <f>. This will modify the conjugated momentum field, but it will 
only contribute a term proportional to 5^i. Again, such a coupling respects the unitarity 
bound. 

With respect to stability, the coupling ( [2.20| ) might alter the dynamics of the Goldstone- 
Stiickelberg field substantially, as it represents a derivative coupling to the matter degrees 
of freedom. However, the coefficient of this derivative coupling is a vector field constructed 
solely from background quantities. Due to the isometries of the Friedman geometry, this 
background vector has to be of the form 



dT 



(1,0,0,0). (2.22) 



Therefore, no spatial derivative enters the coupling ( |2.21| ). Thus the modes experience an 
additional source proportional to dt5^i, and as long as the matter sector takes good care 
of itself this source cannot alter the stability bound, in particular at large proper momenta. 



3. Discussion 



Let us consider the concrete expressions for the coefficients A and B and discuss the possible 
implications of the corresponding bounds. 

3.1 De Sitter spacetime 

For a de Sitter spacetime, A = 6m 2 (m 2 — H 2 ) = —B. The absence of negative norm 
states requires the unitarity bound m 2 > H 2 to hold, which is just the well-known Higuchi 
bound (in our conventions). Provided this bound has been satisfied, classical stability is 
established automatically, as A = — B. In fact, de Sitter geometry is unique within the 
class of Friedman spacetimes where classical and quantum stability requirements coincide, 
and where stability is solely expressed in terms of two model parameters, i.e. the graviton 
mass and the constant de Sitter curvature scale. Since linearisation is permissible, the 
stability bound, m 2 > H 2 , should be an important consistency condition for a nonlinear 
completion, as well. 
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3.2 Generic Friedman spacetimes 

For a generic Friedman spacetime, the absence of negative norm states demands m 2 > 
H 2 + H, while classical stability requires m 2 > H 2 + H/3. We were able to confirm these 
stability bounds by a full-fledged perturbation analysis of massive cosmological gravitons 
involving all degrees of freedom, as well as a complete set of couplings. (See ||].) The same 
bounds were also derived for the special case of scalar field matter in ||12|| . The findings 
of p3| , however, do not coincide with ours, due to the unconventional matter Lagrangian 
that has been used in their paper. 

In the case of an expanding Universe the unitarity bound will always be satisfied 
during radiation (m 2 > —H 2 ) or matter domination (m 2 > —H 2 /2). Thus, the absence 
of negative norm states is guaranteed by the isotropic expansion. This does, however, 
not imply that classical stability is unchallenged during the Universe's expansion, since 
the classical stability bounds m 2 > H 2 /3 (during radiation domination) and m 2 > H 2 /2 
(during matter domination) will eventually be violated when the modes are being evolved 
backwards in cosmological time. 

The situation becomes more interesting for a universe filled with a mixture of matter 
and a cosmological constant. In this case, there is no unitarity issue at early times before 
the transition from matter to cosmological constant domination. However, after the matter- 
cosmological constant transition H 2 + H will become positive and eventually constant as 
the Universe evolves towards its de Sitter fate. If the cosmological graviton's mass is larger 
than the asymptotic value for the Hubble parameter set by the de Sitter curvature scale, the 
theory always respects unitarity. If, however, the mass parameter is chosen to be smaller 
than the cosmological constant, it is guaranteed that the 4> modes first enter the epoch of 
classical instability, before (evolving forward in time) they hit the then would-be unitarity 
bound (since H < 0). As a consequence, massive cosmological gravitons are protected 
against unitarity violation by the background expansion. We conjecture that this sort of 
self-protection could be robust against any nonlinear completion. 

We are well aware that a non-linear completion may suffer from another inconsistency 
problem, the so-called Boulware-Deser ghost fll4| . This definitely requires further investi- 
gations, but there are already first indications that certain non- linear deformations could 
be consistent M, [H| . 

Similar remarks apply to more general source terms. For an expanding universe, 
H < 0, so that the massive cosmological gravitons always first enter the strong coupling 
regime, rendering the would-be unitarity bound fictitious and opening an exciting window 
of opportunity towards a self-protection mechanism. Finally, let us state the situation is 
reversed in a contracting universe, where H > 0. Here, negative norm eft states will show up 
in the weak coupling regime. There is no obvious self-protection mechanism in this case. 
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